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ON  SLOW  VISCO-PLASTIC  FLOW* 
by  William  Prager,  Brown  University 

1.  Introduction^  The  incompressible  visco-plastic  material 

considered  in  the  following  is  a Bingham  solid  [l],  The  mechan- 

ical behavior  of  this  solid  is  most  readily  visualized  in  the 
following  manner.  Let  a Newtonian  viscous  liquid,  a perfectly 
plastic  Mlses  solid  [2],  and  a visco-plastic  Bingham  solid  be  sub- 
jected to  the  same  velocity  strain.  The  stress  in  the  Bingham 
solid  is  then  obtained  by  adding  the  stresses  in  the  Newtonian 
liquid  and  the  Mises  solid. 

Whereas  specific  boundary  value  problems  concerning  the  slow 
flow  of  a Bingham  solid  have  been  discussed  in  the  literature  (see, 
for  instance,  [3]),  few  general  results  concerning  this  type  of 
visco-plastic  flow  seem  as  yet  to  be  available.  In  the  present 
paper  a uniqueness  theorem  and  two  extremum  principles  are  estab- 
lished, 

2.  Notation.  Latin  subscripts  take  the  range  1,  2,  3,  and 
the  summation  convention  operates  on  repeated  subscriptso  The  co- 
ordinates are  rectangular  and  Cartesian,  and  differentiation 
with  respect  to  xj_  is  indicated  by  a comma  subscript  followed  by 
the  subscript  i (f}i  = df/dx^). 

The  following  boundary  value  problem  will  be  considered,  The 
incompressible  Bingham  solid  under  consideration  occupies  the  three- 
dimensional  region  V which  is  bounded  by  the  surface  S.  The  body 
f#rce  (per  unit  of  volume)  is  given  throughout  V,  The  surface 

* The  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Research  (Con- 
tract N7onr-35801). 

**  Numbers  in  square  brackets  refer  to  the  bibliography  at  the 
end  of  the  paper. 


All-84 


2 


traction  is  prescribed  on  the  portion  ST  of  S,  and  the  velocity 
Vj^  on  the  remainder  Sy«  If  ST  = S,  the  given  surface  tractions 
must,  of  course,  satisfy  the  conditions  of  equilibrium,  J T^dS  + 

J F^V  - 0,  and  if  Sy  » S,  the  given  velocities  must  satisfy  the 
condition  of  incompressibility,  J'  v^n^dS  = 0,  where  ni  is  the  unit 
exterior  normal  of  S.  From  the  data  on  the  surface,  the  velocity 
field  v^(x)  and  the  stress  field  djj(x)  are  to  be  determined 
throughout  V under  the  assumption  that  the  effects  of  inertia  are 
negligible  when  compared  to  the  effects  of  the  viscous  stresses 
and  yield  stresses. 

It  is  convenient  to  write  the  stress  tensor  in  the  form 

dij  = “ P6ij  + sij> 

where  p = - d^/3  is  the  mean  pressure,  and  sjj  the  stress  devia- 
tion. The  velocity  strain  is  defined  by 

elj  = 2 (vi,J  + (2) 

The  incompressibility  of  the  solid  requires  that 

eii  = vi,i  = °-  (3> 

The  following  positive  in  variants  of  velocity  strain  and 
stress  deviation  prove  useful  in  the  analytical  description  of  the 
mechanical  behavior  of  the  solid: 

I = (2tiJ  Ell)1/2,  J = (1  si3s13)1/2.  (4) 

3«  Basic  relations.  The  general  relations  between  the  stress 
deviation  and  the  velocity  strain  in  an  incompressible  Bingham 
solid  were  first  given  by  Hohenemsor  and  Prager  [4],  [5]«  In  the 
present  notation,  these  relations  have  the  form 


0 if  J < k, 

= A , kx 

] (1  - -J)  Sij  if  J > k. 
By  squaring  <5)  and  using  (4),  we  find  that 


(6) 


U 
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" 0 if  J < k, 

J » k if  J > k« 

Substitution  of  (6)  into  the  second  Eq*  (5)  yields 

sii  s 2(4  + x)£iJ  if  1 * °*  (7) 

Of  course,  whenever  the  stress  remains  below  the  yield  limit,  the 
velocity  strain  vanishes,  and  the  stress  cannot  be  expressed  in 
terms  of  the  velocity  strain. 

Since  the  inertia  effects  are  neglected,  the  stresses  must 
satisfy  the  equations  of  equilibrium;  these  are 

J + Fi  55  0 in  v>  (8) 


and 


*±3*3  = Ti 


on  Sr 


(9) 


Introducing  (1)  into  (8)  and  eliminating  p by  cross  differen- 
tiation we  obtain 

sij,jk  " skj , ji  + Fi,k  " Fk,i  = 0 (10) 

as  the  condition  of  equilibrium  for  the  stress  deviation. 

In  the  following,  the  principle  of  virtual  work  will  be  used 

ii 

repeatedly.  Let  vi  be  a continuous  velocity  field  with  piecewise 
continuous  first  derivations  satisfying  (3)  and  0^  an  entirely 
unrelated  continuous  stress  field  with  piecewise  continuous  first 
derivatives  satisfying  (8),  Furthermore,  set  ei‘«  = l(v V ,•  + v'J  .) 
and  Tj_  = The  principle  of  virtual  work  is  then  expressed 

by  the  equation 

• ••  f r it  f ti 

Ti  v-i  dS  + F.  v,  dV.  (1 1 ) 

Xx  vjxj.  v - J 


eiJ  dV  = 


4.  Uniqueness.  A stress  field  cs^,  a velocity  field  vl5  and 

the  associated  field  of  velocity  strain  £, . = 2 (vH  . + v^  4)  will 

1 j 2 1 > J J j 1 

be  said  to  constitute  a solution  of  the  boundary  value  problem  for- 
mulated in  Sec.  2 if  they  satisfy  the  boundary  conditions  and  Eqs. 
(3)5  (5),  and  (8). 
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Let  us  assume  that  the  boundary  value  problem  admits  a second 
solution  j 5 v]_9  ejj*  Since  each  solution  corresponds  to  the  same 
body  forces  F^,  the  principle  of  virtual  work  yields 


Cc 


- o 


- e 


ij 


)dV  = P (T±  - TjKvi  - v|)dS.  (12) 

J 


The  right-hand  side  of  (12)  vanishes  because  the  first  factor  of 
the  integrand  vanishes  on  Sj  and  the  second  factor  on  Sy.  To  prove 
that  our  Boundary  value  problem  defines  a unique  field  of  velocity 
strain,  we  shall  prove  that  the  Integral  on  the  left-hand  side  of 
(12)  is  positive  unless  throughout  V. 

On  account  of  the  incompressibility  of  the  solid,  the  inte- 
grand on  the  left-hand  side  of  (12)  can  be  written  as  A.  = (s^-s^j) 
*G"ij  - e^j).  In  discussing  the  sign  of  this  expression,  we  must 
consider  four  types  of  region. 

(.1)  Where  neither  ej_j  nor  ejj  vanishes,  we  have  by  (7)  and 

(4), 

A = 2p-(  ei  j - eij  - eid)  + JMl  + I')(II’  - 26^6^). (i3) 

The  first  expression  on  the  right-hand  side  of  (13)  is  positive 
unless  = e{j,  and  the  second  expression  is  non-negative  by  the 
Schwarzian  ineouality.  Thus,  A can  vanish  in  this  type  of  region 


only  if  eij  = e 


i 

i r 


and  hence 


(y^_)  Where  ^ 0 but  = 0,  we  have,  according  to  (5) 

(14) 


J > k,  J'  < k 


A - (s1;)  - sh)eij  - ^ CsijsiJ 


" SljSij 


Here,  the  factors  in  front  of  the  bracket  are  positive.  Moreover, 


siJsiJ 


= 2J  > 2JJ» 


by  (14),  and  hence 

sijsij  - = 


1JS1J  >2JJ' 


sijsij 


(16) 


by  the  Schwarz ian  inequality.  Thus,  A cannot  vanish  in  this  type 
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of  region. 

( ill)  It  can  be  shov/n  in  a similar  manner  that  A cannot 
vanish  where  = 0 but  £ 0, 

(fv)  Where  = 0,  we  have  A = 0. 

It  follows  from  the  preceding  discussion  that  the  left-hand 
integral  in  (12)  can  vanish,  as  it  do,  only  if  = ell 

throughout  V.  The  field  of  velocity  strain  is  therefore  seen  to 


be  unique.  Since  this  field  determine*?  the  velocity  field  to  with- 
in a rigid  body  motion,  the  velocit-  field  will  also  be  unique 
except  when  Sy  = 0.  In  this  cate,  cha  velocity  field  is  deter- 
mined only  to  within  a rigid  body  motion.  As  regards  the  stress 
field,  this  is  obviously  not  unique  in  such  regions  where  the 
velocity  strain  vanishes.  Where  the  velocity  strain  does  not  van- 
ish, however,  we  have  Sjj  = The  corresponding  stress  fields 

Ojj  and  therefore  differ  at  most  in  the  mean  pressure.  Since 
each  of  these  stress  fields  must  satisfy  the  equations  of  equili- 
brium (8),  this  difference -in  mean  pressure  must  be  constant 
throughout  V.  The  boundary  conditions  on  S*p , however,  rule  out 
such  a difference  in  mean  pressure.  Thus,  the  stress  field  is 
unique  except  when  = 0.  In  this  case  the  stress  field  is  deter- 
mined only  to  within  a constant  hydrostatic  pressure. 

5.  Minimum  principle  for  velocity  strain.  A continuous 

jl; 

velocity  field  vj_  with  piecewise  continuous  first  derivatives  will 
be  called  kinematically  admissible  if  it  satisfies  the  boundary 
conditions  on  Sy  and  the  condition  of  incompressibility  (3),  Where 
the  corresponding  velocity  strain  e*j  does  not  vanish,  Eq.  (7) 
furnishes  a field  of  stress  deviation  s*j.  As  a rule,  however, 
this  field  will  not  satisfy  the  condition  of  equilibrium  (10). 
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Excluding  the  case  Sy  “ 0,  we  propose  to  compare  the  unique 
solution  of  our  boundary  value  problem  to  a generic  kinematic 
cally  admissible  velocity  field  v*«  Consider  the  expression 

H*  = + 2kX*)dV  - 2 J ^vJdSj  - 2 J Ftv*dV  (17) 

which  is  a functional  of  the  velocity  field  v|.  Since  v*  = on 
Sy,  we  have 


H*  - H = 

- 2 


*2 


i U(I  - I ) + 2k(I*  - I)  ]dV  - 

F1(v!  - v.,)dV, 


(18) 


Tl(vi  “ Vi)dS  - 2 
where  the  surface  integral  is  extended  over  the  entire  surface  S. 
We  shall  prove  that  H*  - H >0  unless  v^  = v^,  thus  establishing 
the  following  ixnimum  principle. 

Theorem  1,  Among  all  kinematically  admissible  velocity 
fields,  the  actual  velocity  field  minimizes  the  expression  (17). 

Using  the  principle  of  virtual  work  and  the  condition  of 
incompressibility,  we  write  Eq.  ( 18 ) in  the  form 
H*  - H = J [>(E*2  - I2)  + 2k( I*  - I)  - 2slj(eJJ  - e^JJdV,  (19) 
where  sij  is  the  (unique)  stress  deviation  associated  with  the 
solution  v^. 

In  discussing  the  sign  of  the  integrand  B in  (19)  we  must 
distinguish  two  “types  of  region. 

(j.)  Where  £ 0,  we  have  by  (7)  and  (4) 

B = 2^(e:j  - + ^(11*  - 2 £ i j e ) . (20) 

The  first  term  cm  the  right-hand  side  of  (20)  is  positive  unless 
% 

eij  = eij,  4ie  S5Gond  term  is  non-negative  by  the  Schwarzian 
Inequality,  Thus,  B is  positive  in  region  (i)  unless  e^j  =■  e*^, 
(ii)  When  = 0,  the  integrand  B of  (19)  reduces  to 

B = pil*2  + 2kl*  - 2SijeJj,  (21) 

where  s^j  is  no  Longer  uniquely  determined  but  must  satisfy  the 
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% 


k 


(22) 


according  to  (5).  If  = 0,  the  expression  (21)  vanishes.  If, 
on  the  other  hand,  t 0,  the  associated  stress  deviation  s-jj 
is  readily  found  for  (7)  and  we  have 

J*  > k (23) 

in  accordance  with  (5).  The  expression  (21)  can  then  be  written 
as  follows: 

B = + k - -Uil),  (24) 

4 J* 

Nov/,  < 2 JJ* , by  the  Schwarz! an  inequality.  Equation  (24) 

therefore  yields 


B > ii-=J£(  J*  + k - 2J)  > 0 
4 


(25) 


by  (22)  and  (23). 

It  follows  from  the  preceding  discussion  that  H*  > H unless 
eiJ  = eiJ  throughout  V8  In  view  of  the  boundary  conditions  on 
Sy  this  means  that  H*  > H unless  v*  = throughout  V.  This 
establishes  the  minimum  principle  of  Theorem  1,  In  the  special 
case  of  plane  flow  and  under  the  assumption  that  the  boundary  con- 
ditions involve  only  surface  tractions  and  are,  moreover,  such  as 
to  cause  visco-plastic  deformation  throughout  Vs  an  equivalent 
variational  principle  has  been  given  by  Ilyushin  [6], 

8.  Maximum  principle  for  stress.  A continuous  stress  field 
°Ij  with  piecewise  continuous  first  derivatives  will  be  called 
statically  admissible  if  it  satisfies  the  boundary  condition  on 

Sp  and  the  equations  of  equilibrium  (8),  By  means  of  (5)  a.  field 

* * 
of  velocity  strain  is  associated  with  the  stress  field 

As  a rule,  however,  this  field  of  velocity  strain  does  not  satisfy 

the  compatibility  conditions  and  hence  cannot  be  derived  from  a 

velocity  field. 

Excluding  the  case  S«p  = 0,  we  propose  to  compare  a solution 
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d . . of  our  boundary  value  problem  to  a generic  statically  adraissi- 

**■  tJ 

ble  stress  field  d.jj,  It  should  be  kept  in  mind  for  the  following 
that  two  solution  stress  fields  can  differ  only  in  the  common 
rigid  region. 

Consider  the  expression 

K*  = 2 J T*vidSv  - jl.  J [|j*  - k!  + J*  - kfdV  (26) 

which  is  a functional  of  the  stress  field  Since  on 

S^,  we  have 

K - K*  = 2 J (Ti  - T*) vidS 

" T^j  |5|j  - ki  + [ | J*  - k | + J*  - kfjdV,  (27) 


where  the  first  integral  is  extended  over  the  entire  surface  S. 

We  shall  prove  that  K > K*  unless  is  also  a solution  stress 
field  in  which  case  K = K*.  The  following  maximum  principle  will 
thus  be  established. 

Theorem  2.  Among  all  statically  admissible  stress  fields, 
the  solution  stress  fields  maximize  the  expression  (26). 

Using  the  principle  of  virtual  work  and  the  condition  of  in- 
compressibility, we  write  (27)  in  the  form 

K - K*  =1  - SiP  - i[lj  - kl  + J - kf 

^ 2 "i 

+ J[|J*  - k|  + J*  - k]  j dV,  (28) 
where  is  the  (unique)  velocity  strain  associated  with  the 
solution 

In  discussing  the  sign  of  the  integrand  C in  (28)  we  must 
distinguish  four  types  of  region, 

(i)  Where  J > k and  J*  > k,  we  have  by  (5)  and  (4) 

C = 2 (siJ  - - *ip  - ^(2JJ*  - t2?) 

The  expression  inside  the  parenthesis  in  the  second  term  on  the 
right-hand  side  of  (29)  is  positive  unless  = cs-jj  in  which 
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case  it  vanishes.  In  this  case,  C > 0 if  c / I and  C ~ 0 if  c = 1 

In  all  other  cases  C Is  shown  to  be  positive  by  the  following  argu 


oient.  Since  J > k in  region  (j,),  we  have 

C >^(sirs|j)(sirs^)-(2JJ^»si;)s*;j)=(J»J«)2  >0.  (30) 

(11)  Where  J > k but  J*  < k we  cannot  have  s*  ^ = S-jj.  The 

integrand  C of  (28)  can  be  transformed  as  follows  by  means  of  (5) 

and  (4):  , o * p 

C = (1  - |)(2J  - si;JSij)  «*  (J  - k r 

= (J2  + ^ _ si-jsjj).  (3D 

The  factor  in  front  of  the  parenthesis  on  the  right-hand  side  of 
(31)  is  positive  in  region  (jjj.  Moreover,  since  J > k > J* , 

J2+kJ-sijS*j  > <I2+J*2-sijS^j=  iCs^-s^Ks^-s^j)  > 0.  (32) 

Thus,  C > 0 in  region  (11)* 

(ill)  Where  J < k but  J*  > k,  we  have  = 0,  The  integrand 
C in  (28)  thus  assumes  the  form 

C = (J*  - k)2  > 0.  (33) 

(iv)  Where  J <k  and  J*  < k,  the  integrand  in  (28)  vanishes. 
It  follows  from  the  preceding  discussion  that  K > K*  unless 
s^j  = Sjj  in  all  of  V but  the  rigid  region  (Iv).  This  means  that 
K > K*  unless  s*j  Is  the  stress  deviation  field  associated  with 
a solution.  Theorem  2 Is  thus  established. 

7*  Related  extremum  principles.  The  Newtonian  viscous  fluid 
and  the  perfectly  plastic  Mises  solid  are  limiting  cases  of  the 
Bingham  solid  considered  here.  Indeed,  with  k = 0,  Eq.  (5)  charac- 
terizes a Newtonian  liquid  and,  with  p - 0,  Eq.  (7)  characterizes 
a Mises  solid.  Accordingly,  our  extremum  principles  should  contain 
the  extremum  principles  for  the  Newtonian  liquid  and  the  Mises 
solid  as  special  cases. 

With  k ~ 0,  the  functional  (17)  reduces  to  an  expression  which 
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is  readily  recognised  as  the  excess  of  th#  rate  of  dissipation  of 
energy  in  the  Newtonian  liquid  over  double  the  rate  at  which  work 
is  being  done  by  the  body  forces  and  the  surface  tractions  on  ST. 
That  thij  quantity  is  minimised  by  the  actual  velocity  field  was 
observed  by  Helmholtz  [ 7] • 

With  p ~ 0,  on  the  other  hand,  the  functional  (1?)  reduces  to 
an  expression  which  is  easily  identified  as  double  the  excess  of 
the  rate  of  dissipation  of  energy  in  the  Miscs  solid  over  the  rate 
at  which  work  is  being  done  by  the  body  forces  and  the  surface 
tractions  on  S^.  That  this  quantity  is  minimized  by  the  actual 
velocity  field  was  recognized  by  Harkov  [8]. 

In  view  of  the  fact  that  these  special  cases  of  our  minimum 
principle  can  be  stated  in  terms  of  energy  dissipation,  it  is 
worth  noting  that  the  first  integral  on  tho  right-hand  side  of 
(17)  does  not  represent  the  rate  of  dissipation  of  energy  in  the 
Bingham  solid.  Inueed,  Eq,  (7)  shows  the  rate  of  energy  dissipa- 
tion to  be  r f o 

\ 3i  j ei;j  dV  = J (Hi  + fcl  ) dV.  (34) 

As  Ilyushin  [6]  has  remarked  with  reference  to  the  special  case 
considered  by  him,  the  first  integral  on  the  right-hand  side  of 
(17)  can  be  interpreted  as  the  sum  of  the  "power  of  deformation 
and  the  power  of  the  internal  forces  of  plastic  resistance".  Alter- 
natively, this  integral  can  be  considered  as  the  sum  of  the  rate 
of  viscous  dissipation  of  energy  and  double  the  rate  of  plastic 
dissipation  of  energy. 

With  k = 0,  the  functional  (26)  reduces  to  an  expression  which 
represents  the  excess  of  double  the  rate  at  which  work  is  being 
done  by  the  surface  tractions  on  Sy  over  the  rate  of  dissipation  of 
energy  in  the  Newtonian  liquid,  the  latter  rate  being  expressed  in 
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The  principle  that  this  expression  is  maximized 
by  the  actual  stress  field  seems  to  be  new. 

Finally,  the  characteristic  equation  for  a Mises  solid  in 
plastic  flow  is  obtained  from  (5)  by  letting  p and  J - k tend  to 
zero  in  such  a manner  that  2pJ/(J~k)  tends  towards  a positive 
factor  of  proportionality  between  s.y  and  If  this  limiting 

process  is  applied  to  the  expression  (26),  this  is  seen  to  reduce 
to  double  the  rate  at  which  work  is  being  done  by  the  surface  trac- 
tions on  Sy.  The  fact  that  this  rate  of  work  is  maximized  by  the 
actual  surface  tractions  was  recognized  by  Hill,  first  [9]  in  the 
special  case  where  the  boundary  conditions  are  such  as  to  ensure 
plastic  flow  throughout  V,  and  later  [ 10]  in  the  general  case. 
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